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Methods for calculating the bulk modulus of glasses with homogeneous and heterogeneous structures are pro¬ 
posed within a model that takes account of the different degrees of interaction between structural elements 
comprising two-component glasses. Lead-silicate and lead-borate glasses are studied as examples of the appli¬ 
cation of these computational methods. 
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The existing methods of determining the elastic proper¬ 
ties of solids by measuring the ultrasonic wave velocities 
(longitudinal and transverse) presuppose uniformity of the 
medium in which the waves propagate, i.e., all parameters 
which are measured are constant: the density of the wave- 
propagation medium, wave velocity, and elastic properties. 
For a heterogeneous structure the measured value of the 
wave velocity is a result of averaging, whose algorithm is un¬ 
known, as is whether or not the characteristics are averaged 
according to the same algorithm as the velocities. Thus, the 
application of this method to glass with heterogeneous struc¬ 
ture is fraught with possible errors associated with incorrect 
application of the measured parameters — density, sound ve¬ 
locity — if the sizes 5 of the nonuniformities are comparable 
to the wave length X and completely justified if X » 5. 

In the present work a method is proposed for calculating 
the elastic properties of glasses with a nonuniform structure 
on the basis of data on a property of glasses whose composi¬ 
tion corresponds to the limits of this range, i.e., it corre¬ 
sponds to uniform structures. 

We shall study two-component glasses with a nonuni¬ 
form structure and different contents A and B of the glass for¬ 
mer and modifier, respectively, in glass compositions of the 
type xB - (I —x)A, where x is the molar fraction of the modi¬ 
fier with 0 < x < x x and Xj is the limit of the composition 
range for glass with nonuniform structure. As proposed in 
[1,2], in this range the phases with constant composition be¬ 


have as structural elements (SE) whose properties are inde¬ 
pendent of the composition. The glass former A is denoted as 
SE1, and the other structural element SE2 is the phase whose 
composition is x 1 B - (1 - x 1 )A. 

Since the nonuniformity of these glasses is due to the 
particulars of the interaction between the structural elements 
(the interaction between SE of the same kind is much stron¬ 
ger than that between SE of different kinds), it can be shown 
[2] that the glass density p can be represented as 


P = 


Pi + 


(P2 ~Pl )* 

k(x\ —x) + x 


( 1 ) 


where k = V l /V 2 is the ratio of the molar volumes of the cor¬ 
responding structural elements, while p, and p, are the den¬ 
sities corresponding to the structural elements SE1 and SE2. 
The expression (1) can be put into the form [1] 

p = ax 2 + bx + c, (2) 


where a, b, and c contain the same parameters of the SE but 
here we shall employ the initial expression (1). 

It follows from the dependence of the relative change of 
the density on the pressure dp under hydrostatic compression 


dp dp 
P ~ K 


( 3 ) 
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p, g/cm 3 



Fig. 1. Glass density versus the modifier PbO content: 
1 ) PbO-Si0 2 , L. I. Denikina’s data presented in [3]; 2 ) PbO-B^C^: 
x<0.2, data of R. R. Shaw and D. R. Uhlman presented in [5]; 
x > 0.2, M. S. Heinrich’s data presented in [4], 


show using the expression (1) that the density change accom¬ 
panying a pressure change is 

d p = dp + (dp 2~ dp i^ (P2 -Pi )(*i ~ x )x (4) 
k(x l — x) + x \k(x\ — x) + x] 2 


Introducing the variable p = p /K, we obtain 


dp = ydP; 

x _ p-p t 


k(x l -x)+x p 2 — p i 
A-Oi -x) p 2 -p 


£(x, -x)+x p 2 -p[ 
^dp 2 dpi 

, P 2 Pi 


d/c = k 


( 5 ) 


For two-component glasses with a uniform structure it is 
convenient to use the expression (2) in which 

c = Pf 

6 = (P 2 -Pi) + Pi 2 ; 
a p 12 . 

Here, the initial components can play the role of the SE. 
Since the properties of each component change when the 
components interact with one another, in the linear approxi¬ 
mation of the dependence of the properties of one component 
on the content of the other component the parabolic depend¬ 
ence (2) is obtained for the density. In this case the structural 
elements SE1 and SE2 interact, forming a phase with vari¬ 
able composition, and p 12 is an indicator of the interaction 
between the structural elements SE1 and SE2. Assuming the 
compositions of the SE coincide with the components, x is 
the SE2 fraction in the solid solution SE1 + SE2. 

In contrast to glasses with a nonuniform structure, the 
concept of local density is introduced for the glasses being 
studied here: p, — when a SE is surrounded by SE which are 
similar to it and p' — when a SE is surrounded by elements 
of a different kind. As a result 

Pi 2 = Pi -Pi+P' 2 -P 2 = A Pi +A P 2 - (?) 

Having tied the density change dp to the pressure incre¬ 
ment d/>, we obtain the function y(x ) 

Y = a y x 2 + b y x + c y , (8) 

where 

_ f dp 2 dp l \ 

V P 2 Pi ) , 
by = Y 2 -Yi -a T ; 

Cy =Yi, 


and using the relations (3) and (5) we write instead of (4) 


Y = Yt 


+ P~Pi 
P 2 “Pi 


Y 2 “Yi +(P 2 


-P) 


' 1 

1 Y 


k 2 J_ 


( 6 ) 


where , y 2 and K i , K 2 are the values of the corresponding 
quantities at the limits of the composition range of nonuni¬ 
form structures. 

Since the compositions at these limits are uniform and 
correspond to the structural elements SE 1 at x = 0 and SE2 at 
x = x 1; no errors arise on averaging over the sound propaga¬ 
tion velocities in a nonuniform medium. 

Thus, knowing the characteristics of the nonuniform 
compositions at the range limits they can be calculated for 
any composition inside the range. 


and Y|, Y 2 an d K\, K 2 are the values of y and K for the corre¬ 
sponding structural elements of the glasses. In the case at 
hand such SE are the oxides of silicon and lead for the 
lead-silicate glasses and the oxides of boron and lead for 
lead-borate glasses. 

Thus, in the composition range of uniform structure 
(xj < x < 1) the dependence y(x), just as p(x), is parabolic. 

As an example of the application of the method, we shall 
consider lead-silicate and lead-borate glasses. It follows from 
the form of the composition dependence of the density of 
these glasses (Fig. 1 — sections of the computed curves with 
positive curvature) that in some content range for the modi¬ 
fier PbO (0 < x < x 1 ) the structure of these glasses is hetero¬ 
geneous; for silicate glasses x { = 0.5 molar fractions while 
for borate glasses Xj = 0.2 molar fractions. The data for the 
modulus K and the computed curves for these glasses are 
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presented in Fig. 2. The values of the computed parameters 
are indicated in Table 1. 

Handbook data for quartz, lead-silicate, and lead-borate 
glasses were used for the calculations [3,4]. Preference was 
given to works where the density and elastic characteristics 
were determined for one and the same samples. Comparing 
the figures shows that the spread of the experimental points 
in determining the density is much smaller than in determin¬ 
ing the bulk modulus. 

The parameters of the computed functions p(x) and 
K (x) in the composition range x < x t were determined using 
the experimental values of P] = p(0), K x = K( 0), p 2 = p(x 3 ), 
K 2 =K(x 1 ) and the computed values of k = V l /V 2 , where 
k= 1.26 for lead-borate glass and k = 1.17 for lead-silicate 
glass. Next, the relation (6) was used to calculate Y = Y(p) 
andK(p) = p/y. 

The function y(x ) following the parabolic law (8) in the 
composition range x > x l was determined from the experi¬ 
mental data on K (x) for lead-borate glasses and the com¬ 
puted values of p(x). This made it possible extrapolate y(x ) 
to the entire composition range x l < x < 1 and then compute 
the function K(x) = p(x)/y(x). 

For lead-borate glasses, the individual computed curves 
in the composition range with homogeneous and heteroge¬ 
neous structure are in quite good agreement with the experi¬ 
mental points but do not match with one another near x = Xj. 
This can be explained by the fact that in the composition 
range near a boundary separating glass with homogeneous 
and heterogeneous structure there is also no linking between 
the experimental data (they were obtained by different au¬ 
thors and by different methods: R. Shaw and D. Uhlman used 
echo pulses, A. Osaka and colleagues used the pulse super¬ 
position method), while the computed curves are constructed 
directly from the experimental data. 

For lead-silicate glasses, the data on the elastic properties 
in the composition range for homogeneous structure are in¬ 
adequate to establish from them the dependence y(x), while 
for glasses with nonuniform structure the results of B. I. Ki- 
sin and 1. N. Korolenko from S. V. Nemilov’s work, which 
are presented in [4], are characterized by a small spread. 
These data were used to determine the 
parameters of the curve for x<Xj. The 
parameters of the curve y(x ) for the com¬ 
position range x > Xj were calculated ac¬ 
cording to the values of y(0) and y(0.5) 
for lead-silicate and y(l) for lead-borate 
glasses, i.e., the calculation of these pa¬ 
rameters is in no way related with exist¬ 
ing experimental data (only two points) 
and is based on the suppositions that the 
properties of SE2 for x = 1 are the same 
for lead-borate and lead-silicate glasses, 

SE1 is common to both glasses with ho¬ 
mogeneous and heterogeneous structure, 
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Fig. 2. Bulk modulus K versus the content x of the modifier PbO: 

1) PbO-SiCL: la) [4] (S. V. Nemilov cites measurements performed 
by B. I. Kisin and I. N. Korolenko); lb ) measurements performed 
R. R. Shaw and D. R. Uhlman, presented in [3]; lc) curve computed 
according to the values of K at the points x = 0, x = 0.5, and x = 1; 

2) PbO-Bo0 3 : 2a) x<0.2 according to data of R. R. Shaw and 
D. R. Uhlman, presented in [3]; 2b) x> 0.2, according to data ob¬ 
tained by A. Osaka, N. Soga, and M. Kunuji, presented in [4], 


while the elastic properties of glasses from both ranges are 
the same at x = 0.5. 

Figure 2 displays the computed curves for two types of 
glasses with different glass formers and the same modifying 
agent. 

The calculation of the curve lc relied on the data for pure 
silica and lead oxide as well as on the extrapolation of the 
computed curve la to the composition x = 0.5. The data for 
lead oxide were obtained by extrapolating the curve 2b to 
x — 1. Figure 2 also displays the experimental points, ob¬ 
tained by R. Shaw and D. Uhlman and published in the hand¬ 
book [3], which are not taken into account at all and only 
demonstrate the spread of the data obtained by different au¬ 
thors (curve la and points lb). 

Thus, the satisfactory agreement between the computed 
curves and the experimental points in the composition range 
for nonuniform structure indicates that the average size of 
the nonuniformities 5 <c X ~ 100 pm. On the other hand, for 


TABLE 1. Values of the Computed Parameters 


Glass 

Jtj, molar 
fractions 

Pi. 

g/cm 3 

P2. 

g/cm 3 

Pl2’ 

g/cm 3 

K x x 10 ~ 10 , 
Pa 

K 2 x 10 - 10 , 
Pa 

a y x 10 10 , 
g/(cm 3 • Pa) 

PbO-SiO, 

0.5 

2.2 

6.05 

- 

3.28 

3.61 

- 


0.5 

2.2 

9.01 

2.43 

3.28 

1.96 

3.7 

PbO-B^ 

0.2 

1.84 

3.35 

- 

1.33 

4.15 

- 


0.2 

1.84 

8.95 

0.93 

1.33 

1.98 

6.66 


Note. The numerator and denominator refer to glasses with heterogeneous and homogeneous 
structure structures, respectively. However, for the glasses PbO-Si0 2 , p 2 = 6.05 g/cm 3 and 
K 2 = 3.61 x 10 10 Pa are the density and modulus K of the glass at the heterogeneous structure 
limit, i.e., x = x 1 = 0.5 molar fractions. 
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this spread of the measurements of the elastic properties ob¬ 
tained by different authors by different methods, the method 
proposed here cannot be used as a theoretical method of cal¬ 
culating these characteristics, because it depends on the re¬ 
sults of the measurements. But the theoretical curves y(p) ob¬ 
tained for homogeneous and y(x) for heterogeneous glass 
composition ranges give a rational method of “smoothing” 
the existing experimental results. 
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